Phase delay time and superluminal propagation in barrier tunneling 
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In this work we study the behaviour of Wigner phase delay time for tunneling in the reflection 
mode. Our system consists of a circular loop connected to a single wire of semi-infinite length in 
the presence of Aharonov-Bohm flux. We calculate the analytical expression for the saturated delay 
time. This saturated delay time is independent of Aharonov- Bohm flux and the width of the opaque 
barrier thereby generalizing the Hartman effect. This effect implies superluminal group velocities 
as a consequence. We also briefly discuss the concept called "space collapse or space destroyer". 
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I. INTRODUCTION 

Quantum tunneling is one of the most important phe- 
nomenon with wide range of applications in modern tech- 
nology. In spite of its remarkable success in various ap- 
plications, there still remains the fundamental question, 
namely, how much time does a particle take to traverse 
the barrier? ( so called tunneling time problem). This 
problem has been approached from many different view- 
points, but still there is a lack of consensus about the 
existence of simple expression for this time PI0,E]- This 
is due to the fact that there is no Hermitian operator 
associated with this time. The various time scales pro- 
posed in the literature are based on different operational 
definitions and physical interpretations. They may repre- 
sent different complimentary aspect of tunneling process. 
Some of them include dwell time, Larmor time, complex 
times, Buttiker- Landauer times, sojourn times Q, time 
scale based on Bohm's view phase delay time, etc. 

In our present work we concentrate mainly on the con- 
cept of Wigner phase delay time. This time is usually 
taken as the difference between the time at which the 
peak of the transmitted packet leaves the barrier and the 
time at which the peak of the incident Gaussian wave 
packet arrives at the barrier. Within the stationary phase 
approximation the phase time can be calculated from the 
energy derivative of the 'phase shift' in the transmitted 
or reflected amplitudes. This phase delay time is also 
related to various physical quantities such as partial den- 
sity of states |6|. In the case of a quantum tunneling 
it has been shown that in the opaque barrier limit the 
phase delay time does not depend on the barrier width. 
This phenomenon is called as the Hartman effect. This 
implies that for sufficiently large barriers the effective ve- 
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locity of the particle can become arbitrarily large, even 
larger than the speed of light in the vacuum (superlumi- 
nal). This may be further regarded as the evidence of 
the fact that quantum systems seem to behave as non- 
local. Experiments have clearly demonstrated @, 0, 
that "tunneling photons" travel with superluminal group 
velocities. Their measured tunneling time is practically 
obtained by comparing the two peaks of the incident and 
transmitted wave packets. It is important to note that 
these observations do not violate 'Einstein causality', i.e., 
the signal velocity or the information transfer velocity is 
always bounded by the velocity of light. Theoretically 
Hartman effect has been generalized to different cases 
including double barriers, various geometric structures 
and in the presence of Aharonov-Bohm flux. Mystery 
surrounding this effect has also been addressed recently 
by Winful. He argues that the short time delay observed 
is due to energy storage and release and has nothing to do 
with propagation. Since the stored energy or the proba- 
bility density in the evanescent field decreases exponen- 
tially with the distance in the barrier, after certain decay 
distance it does not matter how much more length the 
barrier has. For details we refer to [Toj . 

In our present work we study the above mentioned ef- 
fect in a simple quantum system consisting of a circular 
loop connected to a single wire of semi-infinite length as 
shown in Fig^ I n addition we impose Aharonov-Bohm 
flux through the loop. We also focus on the following 
situation. There is a potential barrier (or barriers) inside 
the loop, while the potential in the connecting lead is set 
to zero. The incident energy of the free propagating elec- 
tron E in the semi-infinite wire is less than the barrier 
potential height V. The impinging electron on the loop in 
the sub barrier region travels as an evanescent wave in the 
loop before being reflected. We analyze the phase time of 
the reflected wave. We show that this phase time in the 
opaque barrier regime becomes independent of the length 
of the circumference of the ring and the magnitude of the 
Aharonov-Bohm flux, thereby generalizing the Hartman 
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FIG. 1: Schematic diagram of a circular loop connected to 
semi infinite lead. 



effect beyond one dimension and in the presence of mag- 
netic flux. We have also extended this effect by including 
an additional potential well in the circular ring. Interest- 
ingly the saturated delay time becomes independent of 
the length of the potential well (in the large length limit) 
in the off resonant case. The potential well can support 
many resonant states. This result is regarded as a "space 
collapse or space destroyer" ^lj . Even though in the po- 
tential well inside the loop electronic wave can travel as a 
free propagating mode (and not as a evanescent mode), 
surprisingly the saturated delay time is independent of 
the length of the well (as if it does not count). 

In the sections to follow we first give the description 
of the theoretical treatment. For this we make use of 
the quantum waveguide approach well known for the 
electronic transport problems. The classical relativistic 
Helmholtz equation for the electromagnetic wave propa- 
gation is formally identical to non relativistic stationary 
Schrodinger equation. Hence conclusions obtained using 
quantum tunneling theory are equally valid in the case 
of electromagnetic phenomenon ( photon tunneling ). In 
the later sections we will discuss our results and conclu- 
sions. 



II. THEORETICAL TREATMENT 



We approach this scattering problem using the quan- 
tum wave guide theory [T2I Il3j. In the stationary case 
the incoming particles are represented by a plane wave 
of unit amplitude. The effective mass of the prop- 



agating particle is m and the energy is E = h k /2m 
where k is the wave vector corresponding to the free par- 
ticle. The wave function in different regions (which are 



solutions of the Schrodinger equation) in the absence of 
magnetic flux are given below 

V, (x ) = e lkx ° + Re- lkxa ( in region I ) (1) 

M x i) = A ie lqiXl + B ie -' lqiXl (in region II ) (2) 

MM = A 2 e lkX2 +B 2 e- lkX2 ( in region III) (3) 

MM = A 3 e' q3 X3 + B 3 e~ lq3 X3 ( in region IV) (4) 

with k being the wavevector of electrons in the lead 
and in the intermediate free space between two barri- 



ers inside the ring. q\ = y 2m(E — V\)/h and q 3 — 

2m(E — V3 ) / 7z 2 are the wavevector respectively for 
propagating electrons in the barriers of strength V\ and 
V3 inside the ring. The origin of the co-ordinates of xq 
and X\ is assumed to be at J and that for x 2 and x 3 
are at P\ and P 2 respectively. At Pi, x\ = Ibx, at P 2 , 
x 2 = w and at J, x 3 — lb 3 , where lb\ and lb 3 are the 
length of the two barriers separated by a well region of 
length w inside the ring. Total circumference of the ring 
is L = lb\ + Z63 + w. 

We use Griffith's boundary conditions [l4] | 



and 



MO) = MO) = Mih) 



dMM , dipi(xi) dip 3 ( x 3) 



dx 



dx\ 



dx 3 



0. 



(5) 



(6) 



at the junction J. All the derivatives are taken either out- 
ward or inward from the junction [T^ . Similar boundary 
conditions are held at Pi and P 2 i.e. 



at Pt 



dx\ 
at P 2 

dip 2 ( x 2) 



+ 



Mlh) 
^2(^2) 



dxi 



dx 2 
^3(^3) 



Pi 



dx 3 



Pi 



Mo) 
= 0, 
^3(0) 
= 0. 



(7) 
(8) 
(9) 
(10) 



We choose a gauge for the vector potential in which 
the magnetic field appears only in the boundar y co ndi- 
tions rather than explicitly in the Hamiltonian |l2l Il5j . 
Thus the electrons propagating clockwise and anticlock- 
wise will pick up opposite phases. The electrons prop- 
agating in the clockwise direction from J will pick up 
phases i a± at P±, i (a± + a 2 ) at P 2 and i (ai + a 2 + a 3 ) 
at J after traversing once along the loop. The total phase 
around the loop is a.\ + a 2 + a 3 = 2tt (fr/tfro, where <j> 
and <j)Q are the magnetic flux and flux quantum, respec- 
tively. Hence from above mentioned boundary conditions 
wc get for tunneling particle 
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1 + R - A± - B\ exp(-i a.\) 


= 0, 


(11) 


A 3 exp(— k 3 lb 3 ) exp(ia 3 ) + B 3 exp(/t 3 lb 3 ) — 1 — R 


= 0, 


(12) 


ik (1 — R) + n\ (Ai — B\ cxp(— i ai)) 






- n 3 A 3 exp(-K 3 lb 3 ) exp(i a 3 ) - n 3 B 3 exp(/t 3 lb 3 ) 


= 0, 


(13) 


A\ exp(— Ki lb\) exp[i a{\ + B\ exp(Ki lb\) 






- A 2 - B 2 cxp(-ia 2 ) 


= 0, 


(14) 


Ki Ai exp ( — K\ lb\) exp(i a-{) — K\ Bi cxp(«; 1 lb{) 






+ ik Ai — ik B 2 exp(— ?a 2 ) 


= 0, 


(15) 


A 2 exp(ifcw) exp(za2) + B 2 e~xp(—ikvu) 






- A 3 ~ B 3 exp(-ia 3 ) 


= 0, 


(16) 


ik A 2 exp(i kw) — ik B 2 exp(— i k w) exp(— ia 2 ) 






+ k 3 A3 — n 3 B 3 exp(— i a 3 ) 


= 0, 


(17) 



with Ki = J2m(Vi-E)/h 2 and k 3 = J 2m(V 3 — E)/h 2 
being the imaginary wave vectors, in presence of rectan- 
gular barriers of strength V\ and V 3 respectively, inside 
the ring. 



III. RESULTS AND DISCUSSION 

Once R is known, the 'reflection phase time' r can be 
calculated from the energy derivative of the phase of the 
reflection amplitude 0,01 as 



= h 



dArg[R] 
dE 



(18) 



where, v = hk/m is the velocity of the free particle. The 
concept of 'phase time' was first introduced by Wigner 
|16| to estimate how long a quantum mechanical wave 
packet is delayed by the scattering obstacle. 

For a ring system with a rectangular barrier of strength 
V\ along its entire circumference we obtain an analytical 
expression for the reflection amplitude as 



R 



—ki (2 cos(a) — exp(fcL)) + i | exp(fcL) 
hi (2 cos(a) — exp(fcL)) + i | exp(fcL) 



(19) 



where a = a\ + a 2 + a 3 . In what follows, let us set 
fi = 1 and 2m = 1. We now proceed to analyze the 
behavior of t as a function of various physical parame- 
ters for different ring systems. We express all the phys- 
ical quantities in dimensionless units i.e. all the barrier 
strengths V n in units of incident energy E (V n = V n /E), 
all the barrier widths lb n in units of inverse wave vector 
k^ 1 (lb n = klb n ), where k — ^/E and the reflection phase 
time r in units of inverse of incident energy E (r = Et) . 
After straight forward algebra in the large length (L) 
limit and in absence of magnetic flux, we obtain an an- 
alytical expression for the saturated phase delay time ( 
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FIG. 2: In absence of magnetic flux (i.e. (f> =0 ), for a ring 
with a barrier of strength V\ throughout its circumference, the 
phase time r is plotted as a function of ring's circumference 
L. The solid, dotted, dashed and dash-dotted curves are for 
Vi = 1.25, 1.5, 3, 5 respectively. Incident energy is set to be 
E = 1. In the inset the saturated value of phase time r s is 
plotted as a function of the barrier's strength for same E. 



using Eq. (|18fl and Eq. (|19fl ), which is given by, 



k ki 



+ 



( 2 + if)' 



(20) 



Here the rectangular barrier has strength V\ and width 
Ib-i = L. 

First we take up a ring system with a single barrier 
along the circumference of the ring. For a tunneling par- 
ticle having energy less than the barrier's strength we 




FIG. 3: For a ring with a barrier of strength Vi lies through- 
out its circumference, the saturated phase time r s is plotted 
as a function of magnetic flux <j>. The solid, dashed and dotted 
curves are for L — 6, 7, 9 respectively. Other system parame- 
ters are Vi = 5, E = 1. 



FIG. 4: In absence of magnetic flux (i.e. <j> =0 ), for a ring 
with two barriers of strength Vi and Vz separated by an in- 
termediate well region, the phase time r is plotted as a func- 
tion of ring's circumference L for different width w of the 
well. The dashed-dotted, solid, dotted and dashed curves are 
for w = 0,1,5,10. Other system parameters are lbs, = 5, 
Vi = V 3 = 2 and E = 1. 



find out the reflection phase time r as a function of bar- 
rier's width L which in turn is the circumference of the 
ring. We see (Fig. 0) that in absence of magnetic flux, r 
evolves as a function of L and asymptotically saturates to 
a value r s which is independent of L thus confirming the 
'Hartman effect'. From Fig- El it is clear that the satura- 
tion value increases with the decreasing barrier-strength. 
In the inset of Fig. [21 we have plotted r s as a function 
of barrier-strength. From this we can see that for elec- 
trons with incident energy close to the barrier-strength 
the value of r s is quite large. 

To see the effect of magnetic flux on 'Hartman ef- 
fect', we consider the same system but in presence of 
Aharonov-Bohm (AB) flux. We find out, for the tun- 
neling particle, the reflection phase time as a function 
of embedded magnetic flux for different lengths L of the 
barrier covering the ring's circumference. We have cho- 
sen the lengths such that in absence of the 'AB-flux', for 
a given system (i.e. for known E and V) the reflection 
phase time r gets saturated in these lengths. From Fig.|3| 
we see that r as a function of </> shows AB-oscillations 
with an average value which is the saturation value t s 
for the same system in absence of AB-flux. Further ob- 
serve that (Fig. |2J r is flux periodic with periodicity 4>o. 
This is consistent with the fact that all the physical prop- 
erties in presence of AB-flux across the ring must be pe- 
riodic function of the flux with a period 4>o EJ El 113 • 
However, we see that as we increase L the magnitude of 
AB-oscillation in r decreases. Consequently in the large 
length limit the visibility vanishes. This clearly estab- 
lishes 'Hartman effect' even in presence of AB-flux. The 
constant value of r thus obtained in the presence of flux 
is identical to r s in the absence of flux (see Fig. J2Jl) in 
the large length regime and its magnitude is given by 
Eq. (|20|) . This result clearly indicates that the reflection 



phase time in the presence of opaque barrier becomes 
not only independent of length of the circumference but 
also is independent of the AB-flux thereby observing the 
'Hartman effect' in the presence of AB-flux. 

Now we consider the ring system with the ring having 
two successive barriers separated by an intermediate free 
space as shown in Fig.^ In absence of magnetic flux, we 
see the effect of 'quantum well' on the reflection phase 
time t. In Fig.0]r is plotted as a function of one of the 
barrier's length (say lb\) while other barrier's length is 
fixed (/&3 = 5) and for few different values of length of the 
well. Here, the fixed value of the barrier's length Z63 is 
chosen in such a way that in absence of the well region the 
reflection phase time reaches saturation at this length. 
From Fig. 2| we see that for all parameter values of well's 
width, the saturation value of reflection phase time r s is 
same and it is equal to what we obtained in absence of 
the well in the ring system. Thus the saturated phase 
time becomes independent of the width of the well (in 
the long length limit) in the off resonant case. This is as 
if the effective velocity of the electron in the well becomes 
infinite or equivalently length of the well does not count 
(space collapse or space destroyer) while traversing the 
ring. 

Finally consider a similar system as that shown in 
Fig. ^ Here we see the effect of resonances, present 
in the ring system with a well, on the saturated reflec- 
tion phase time t s . For the system described above with 
V\ = V3 = 2, Z63 = 5, cf) = we have plotted r s , for the 
electrons with incident energy E — 1, as a function of the 
well's width for different parameter values lb\ in Fig. [3] 
We see that the resonances which have Lorentzian shape 
become sharper and narrower as the width of the bar- 
rier lb\ becomes large. For very large lb\ the resonances 
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=0 ), for 



ring with two barriers of strength Vi and V3 separated by 
an intermediate well region, the saturated phase time r a is 
plotted as a function of ring's circumference L for different 
width lb\ of the barrier. The solid and dashed curves are 
for Ibi = 2.75, 3 respectively. Other systems parameters are 
lb 2 = 5, Vi = V 3 = 2 and E = 1. 

are very hard to detect. It should be noted that as we 
increase the length of the well for fixed E for particular 
barrier lengths incident energy E coincides with reso- 
nances (or resonant states) in the well (which arise due 
to constructive interference due to multiple scatterings 
inside the well). For these values of lengths we observe 
sharp rise in the saturated delay time and its magnitude 



depends on the length of the well. It is worth to mention 
that away form the resonance the value of r s is indepen- 
dent of the length of the well (see Fig0J and depends 
only on the barrier strength. 



IV. CONCLUSIONS 

We have generalized the Hartman effect beyond one 
dimension and in the presence of magnetic flux. This is 
done by studying the phase delay time in a circular ge- 
ometric ring in the reflection mode. We have obtained 
an analytical expression for the saturated delay time. We 
have also shown that in the presence of successive opaque 
barriers separated by the potential well the saturated de- 
lay time becomes independent of the width of the well 
(in the long length limit) in the off resonant case. This 
is as if the effective velocity of the electron between the 
well becomes infinite or equivalently length of the well 
does not count (space collapse or space destroyer). We 
have further extended this effect in quantum mechani- 
cal networks with many arms wherein we have shown 
additional existence of new quantum non-local effects as- 
sociated with the phase delay time 0] . Our reported re- 
sults are amenable to experimental verifications for pho- 
ton tunneling in appropriate electromagnetic geometries. 
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